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ABSTRACT.  The  structural  response  of  a Soviet  TM-46  land  mine  with  two 
viscoelastic  explosive  fills  subjected  to  an  externally  applied  pressure  wave 
has  been  analyzed  with  the  ADINA  finite  element  code.  The  main  charge  consists 
of  5.72  kg  TNT  while  the  booster  charge  in  the  fuze  contains  .04  kg  Tetryl  in 
the  fuze  well.  The  finite  element  model  of  the  mine  uses  the  axisymmetric 
two-dimensional  mesh  configuration  with  a rigid  base  support  boundary  condition. 
Both  implicit  and  explicit  time  integration  schemes  have  been  used  for  this 
analysis. 


The  viscoelastic  explosive  filler  materials  exhibit  marked  nonlinear 
behavior.  It  was  therefore  decided  that' the  tension  cut-off  curve  description 
material  models  were  the  appropriate  models  to  use.  Relationships  between 
the  volume  strains  and  the  bulk  moduli  were  obtained  from  the  Mie-Griineisen 
eguations  of  state.  These  models  include  failure  criteria  which  allow  tension 
cut-off  planes  to  form  in  a direction  normal  to  the  principal  tensile  stress 
whenever  the  strain  initially  exceeds  0.1%  In  tension.  The  materials  for  the 
steel  casing  were  modeled  with  bilinear  stress-strain  curves,  von  Mises  yield 
condition,  and  kinematic  hardening  rule.  Trapped  air  Inside  the  mine  body 
was  modeled  as  an  assembly  of  inviscid  linear  compressible  fluid  elements. 


The  finite  element  model  was  initially  verified  for  mode  shapes  at  a few 
low  order  eigenfrequencies  and  a failure  criterion  for  the  casinq  was  incorpo- 
rated based  on  a comparison  of  the  value  of  the  three-dimensional  second  invari- 
ant of  plastic  strain  with  that  of  the  one-dimensional  value  obtained  from  the 
tensile  tests.  Solution  of  the  problem  In  terms  of  stresses  and  displacements 
out  to  .75  ms  of  real  time  indicates  high  stress  concentration  and  large  dis- 
placements of  the  top  cover  plate  in  the  stepped  region  and  minimal  response  of 
the  interior  structure  until  the  cover  plate  is  in  contact  with  the  intermediate 
partition  of  the  mine  body. 


1.  INTRODUCTION.  This  paper  describes  the  response  of  the  Soviet  TM-46 
Antitank  mine  with  a unique  double  walled  construction  of  the  top  pressure  plate 
designed  to  resist  a transient  blast  load.  The  rationale  for  this  analysis  Is 
the  need  to  develop  a remote,  expeditious  means  of  clearing  a path  through  an 
enemy  mine  field.  A method  of  imparting  a relatively  large  transient  pressure 
and  impulse  to  the  surface  of  the  earth  by  means  of  explosives  Is  under  development 
The  current  study  is  a part  of  a general  Investigation  to  determine  the  extent 
of  structural  damage  sustained  by  the  mine  body  from  a given  level  of  blast  wave 
amplitude  and  shape.  The  principal  kill  mechanism  is  to  be  a serious  distortion 
or  rupture  of  the  mine  casing  rather  than  fuze  Initiation  or  removal  of  pressure 
plate  since  the  damage  mechanisms  could  be  easily  changed  from  a particular  type 
of  mine  to  another  and  a s.urekill  could  not  be  ensured  based  on  a particular  mode 
of  actuation. 
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The  mine  investigated  represents  a typical  Soviet  antitank  mine,  which 
consists  basically  of  a round  thin  metal  body  filled  with  an  explosive. 

The  unique  feature  of  this  mine  is  a double  walled  construction  of  the  top 
wall  resulting  from  joining  of  the  top  pressure  plate  with  the  intermediate 
wall  along  the  circular  edge  through  a crimped  joint  which  tends  to  increase 
the  blast  resistance  behavior  of  the  mine.  The  other  distinctive  feature 
is  the  fuze  mechanism.  However  a variety  of  radically  different  fuzes, 
different  both  in  mechanical  designs  and  method  of  activation  could  be  substi- 
tuted for  the  currently  used  ball-spring  mechanism.  Therefore  the  numerical 
model  in  the  present  study  does  not  include  a model  of  the  fu2e. 

The  paper  has  four  major  areas  as  follows:  (a)  problem  definition,  (b) 

determination  of  material  properties  and  selection  of  failure  criteria  of  the 
casing  as  well  as  numerical  characterization  of  the  viscoelastic  explosive 
fills,  (c)  finite  element  model  description  and  calculations,  and  (d)  dynamic 
response  prediction  of  the  structural  assembly. 

2.  PROBLEM  DEFINITION. 

A.  Antitank  Mine  Description.  The  TM-46  land  mine  has  a cylindrical 
steel  body  with  a primary  fuze  well  in  the  center  of  the  top  and  one  at  the 
bottom,  presumably  for  antilift  or  booby  trapping  purposes.  In  addition,  it 
has  a secondary  fuze  well  in  the  sidewall  underneath  the  carrying  handle.  A 
sectional  drawing  of  the  mine  Is  shown  in  Figure  1.  The  mine  has  a nominal 
diameter  of  29.7  cm,  height  of  7.3  cm  and  weighs  8.7  kg  with  a main  charge 
of  5.7  l,g  TNT. 

The  mine  body  Is  made  of  three  pieces  of  sheet  steel  which  are  joined 
at  the  upper  periphery  by  a 360°  crimp.  The  top  cover  of  the  mine  body  is 
only  0.635  mm  thick  and  has  three  steps.  This  cover  connects  to  a central 
circular  plate  formed  by  spot-welding  of  a thick  plate  to  the  thin  cover  section. 
The  intermediate  wall  is  formed  from  0.94  mm  thick  steel  sheet  to  which  a 
hollow  cylindrical  piece  0.56  mm  thick  Is  attached  to  form  the  centrally  located 
top  fuze  well.  The  fuze  well  contains  a 40  g tetryl  booster  charge  for  fuze 
acti  vation . 

The  lower  part  of  the  mine  body  is  formed  by  a deep  drawing  operation 
which  results  in  very  inhomogenous  material  properties.  The  central  cavity  in 
the  main  body  of  the  mine  is  filled  with  a charge  of  5.7  kg  TNT  explosive.  The 
cavity  between  the  top  and  Intermediate  walls  Is  unfilled.  However  compression 
of  air  In  this  region  can  contribute  to  alteration  of  the  response  behavior  of 
the  mine  and  subsequent  uncrimping  of  the  Joint. 

The  normal  method  of  activation  of  the  fuze  Is  by  means  of  force 
applied  to  the  pressure  cap  depressing  the  fuze  and  releasing  the  striker  to 
strike  the  booster  charge  In  the  fuze  well.  This  activates  the  tetryl  booster 
which  In  turn  detonates  the  primary  TNT  charge.  The  secondary  fuze  well  on 
the  TM-46  mine  gives  It  anti -disturbance  capability. 
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B.  Guidelines  for  the  Numerical  Model,  In  keeping  with  the  objective 
of  identifying  a general  failure  mechanism  independent  of  some  specific  design 
feature,  all  fuzes  and  springs  were  omitted  from  the  finite  element  model  of 
the  mine.  This  was  done  in  accordance  with  the  previously  stated  guideline  of 
not  identifying  failure  of  the  fuze  components.  The  model  considered  for  this 
study  does  not  include  secondary  fuzes  and  filling  holes.  However  the  secondary 
tetryl  booster  charge  is  included  to  facilitate  assessment  of  the  influence  of 
trapped  air  in  the  unfilled  space  below  the  top  wall  on  structural  response  of 
the  mine. 

The  auxiliary  fuze  wells  were  not  considered  in  the  current  investi- 
gation since  they  increase  susceptibility  of  the  mine  to  damage  due  to  stress 
concentration  near  the  junction  between  the  body  and  the  fuze.  Thus,  the 
simplified  model  is  conservative  in  terms  of  blast  load  requirements  for  mine 
deactivation.  Also,  inclusion  of  these  unsymmetrically  loaded  structures  would 
require  the  use  of  a three-dimensional  (3-D)  finite  element  model  resulting  in 
significant  increase  in  computing  time  and  costs.  The  dimples  at  the  base  of 
the  mine  were  eliminated  for  the  same  reasons.  Because  of  these  simplifications 
the  2-D  ax i symmetric  model  was  cost-effective  and  conservative  for  dynamic 
response  evaluation. 

C.  Base  Support  and  Surface  Loading.  During  field  emplacement,  the 
mines  may  be  placed  on  the  surface  and  covered  with  grass  or  other  materials 
for  concealment.  In  other  cases,  the  mine  may  be  shallow  buried.  In  both 
cases,  the  mine  will  experience  transient  pressure  loading  on  the  top  surface 
due  to  detonation  of  a countermine  explosive  in  the  vicinity.  It  is  expected 
that  typical  field  boundary  support  conditions  would  be  bracketed  by  two  extreme 
situations.  In  one  case,  the  mine  could  be  simulated  as  being  buried  in  soil 

up  to  its  top  surface  while  the  base  is  supported  on  nonlinear  springs  as 
described  in  Reference  1.  The  other  support  condition  allows  the  mine  to  be 
supported  on  a rigid  roller  base  which  closely  models  the  experimental  conditions 
described  in  Reference  2.  In  the  current  study  only  the  second  support  condition 
was  simulated  In  the  numerical  model  but  inclusion  of  the  soil  medium  implicitly 
through  nonlinear  spring  node-tie  boundary  elements  or  explicitly  through  dynamic 
property  characterization  of  the  soil  could  be  made  without  significant  change 
in  the  basic  model.  A roller  support  condition  was  used  allowing  lateral,  but 
no  vertical,  motion.  In  this  rigid  support  condition,  the  input  shock  load  Is 
applied  to  the  top  and  sides  of  the  mine. 

For  structural  loading  the  pressure  pulse  used  in  this  paper 
simulated  peak  pressure  and  impulse  measured  from  experiments  conducted  with 
mine  clearance  types  of  explosives  in  Reference  2.  The  peak  pressure  was 
13.8  MPa  and  the  impulse  delivered  was  6.5  kPa-sec.  A decaying  exponential 
function  was  fitted  to  these  parameters  resulting  in  the  following  equation 

P{t)  * 13.76i2117t  (1) 

A curv  c this  function  varying  in  time  is  shown  in  Figure  2. 
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Figure  2.  Shock  Loading  in  Antitank  Pines 

3.  MATERIAL  PROPERTIES  CHARACTERIZATION  AND  FAILURE  CRITERIA.  Material 
properties  were  required  for  the  steel  jackets,  the  explosive  filler  materials 
and  the  trapped  air  before  numerical  simulation  could  be  carried  out.  Mechanical 
properties  were  measured  for  the  steel  jackets  by  employing  uniaxial  tensile 
tests.  The  high  pressure  equation  of  state  data  for  explosives  were  obtained 
from  available  publications.  Failure  criteria  used  for  the  steel  jackets  and  the 
filler  materials  were  similar  to  the  formulations  in  Reference  3. 

A.  Steel  Casing.  The  TM-46  jacket  is  made  of  a low  carbon  soft  magnetic 
steel  equivalent  to  mild  steel.  The  lower  part  of  the  casing  was  deep  drawn, 
but  It  retained  an  equlaxed  grain  microstructure  with  isotropic  properties.  Two 
tensile  speciments  were  cut  from  each  of  the  significant  surfaces  of  the  mine 
body.  Locations  of  these  specimens  are  shown  In  Figure  3(a).  The  specimens 
were  machined  with  a large  radius  on  the  test  section  as  shown  in  Figure  3(b). 

An  extensometer  and  a biaxial  strain  gage  were  attached  at  the  location  of  the 
minimum  width  and  the  specimens  were  tested  in  an  Instron  Testing  Machine  at  a 
relatively  slow  rate  of  strain.  Typical  quasi-static  stress-strain  curves  for 
the  Soviet  mine  body  are  shown  In  Figure  4.  Evidence  of  work  hardening  and 
residual  stress  was  significant  in  the  casing  material  due  to  the  forming 
operation  and  operating  field  conditions. 
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SIDE  WALL 

(a)  LOCATION  OF  SPECIMENS,  TM-46  MINE 


(b)  PREPARATION  OF  SPECIMEN  DIMENSIONS  (cm) 

Figure  3.  Details  of  Tensile  Specimen  Sampling  and  Preparation 
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Bilinear  approximations  to  the  stress-strain  curves  obtained  by  averaging 
the  data  for  the  individual  specimens  are  shown  superimposed  in  Figure  4. 

The  AOINA  (4,  5)  finite  element  code  used  in  this  analysis  has  a bilinear, 
elastic-plastic,  von  Mises  yield  condition,  kinematic  hardening,  axisym- 
metric  2-D  element  for  the  steel  jacket. 

The  criterion  selected  to  predict  failure  of  the  steel  casing 
material  was  described  in  Reference  3 as  the  value  of  the  second  invariant 

p 

of  plastic  deviatoric  strain  at  failure,  ^ e , defined  as 

, , P , P P ,,x 

2f (c  ) ■ **  cij  f ij  (2) 

D 

where  . are  the  plastic  component  of  strains  at  failure.  In  the  uniaxial 
tension  test  where  the  load  is  applied  in  the  axial  Z-direction  we  have, 

>2f('l-D)  ■ 3/4  (4)2  <3' 

B.  Characterization  of  Explosive  fills.  There  are  two  types  of 
explosive  filler  materials  employed  in  the  TM-46  mine,  i.e.,  TNT  as  the  main 
charge  and  tetryl  as  the  fuzewell  booster  charge. 

After  surveying  the  available  material  properties  of  explosives  and 
the  various  2-D  axisymmetric  materials  models  in  the  ADINA  code,  the  curve 
description  material  model  (Section  XII,  pp  17-22,  Reference  4)  was  found 
to  be  the  appropriate  model  to  use.  This  model  requires  tables  of  loading  and 
unloading  bulk  moduli  and  shear  moduli  versus  volumetric  strain. 

A relationship  between  the  volumetric  strain  and  the  bulk  modulus 
obtained  from  the  Mie-Griineisen  equation  of  state  (Reference  1,  6)  for 
shock  propagation  in  solids  is  given  as 

P(u,E)  = An  + Bu2  + Cn3  + rE/V  (4) 

where  A,  B,  and  C are  material  constants  determined  experimentally  from 
Hugoniot-pressure  volume  states  obtained  in  shock  transitions.  The  above 
equation  characterizes  the  P,  V,  and  E state  variables  which  are  attainable 
by  the  solid.  In  this  case  the  solid  Is  either  TNT  or  tetryl  explosive. 

The  adiabatic  loading  bulk  modulus  is  defined  by 


Assuming  that  the  state  of  the  system  is  defined  by  the  variables  s,  E,  and 
V,  we  have  corresponding  to  the  expressions  P = P(V,E)  and  E = E(V,s), 
respectively, 
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dp ' (w)E  d#  4 (i), dE 
d£  ■ (w\  d*  4(§),  ds 

Substituting  dE  from  (7)  into  (6),  we  find 

HD  - IV3p\  a /3p\  /8E\  1 W3P\  /3E\  . 

dP  - [(w)£  * (sf),  (sv)J  dv  4 (5t)#(k)s  ds 

But  from  P = P(V,s)  we  know  that 


dP ' (wl dv  4 (£1 ds 


Comparing  Equations  (8)  and  (9),  we  have  the  following  expression. 


/8P\  _/3P\  . /<*P\  /?E\ 

\®V/S  ~ \ 3V/e  Wv  Ws 


Defining  the  pressure  by  the  relation 


4 - (£1 


and  substituting  this  in  Equation  (9),  we  have 


(lv\  ' (w)E  ■ p (lf)t 


Performing  the  required  differentiation  of  (4),  and  using, 

✓ 3P\  . /DP\  dp 
\3V/e  ' Uu/E  dV 


we  obtain,  after  some  manipulation 


2,  * 


« = Pr  + TE/V  + (1  + u)(A  + 2Bu  + 3CM  ) 

Using  P = P(E,V)  from  (4)  the  relationship  between  bulk  modulus  and 
volumetric  strain  is  given  as 

Kt  = r(r  + 1 ) Ep  + A + A'g  + B ' + C'p^ 


strain  ratio,  u = (V0  - V)/V  and  V0  = 1/P(J. 
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where 

A'  = A(r  + 1)  + 28 
B’  = B(r  + 2)  + 3C 
C'  = C(r  + 3) 

In  order  to  transform  Equation  (15)  to  a form  k = *(V)  the  assumptions  of 
unidimensional  shock  equations  and  conservation  laws  of  mass,  momentum  and 
energy  are  invoked.  The  three  conservation  equations  in  a frame  moving  at 
the  shock  velocity,  Us,  are: 


»’  * povo 

(16) 

»ovoZ 

(17) 

0VE  * S(oV)»2  * PV  . * >,(»„«„)«</  * Vo 

(18) 

where  V and  U are  particles  velocities,  E Is  the  specific  internal  energy 
and  the  quantities  with  zero  subscript  are  undisturbed  values  while  the 
quantities  without  subscripts  are  applicable  behind  the  shock.  In  a 
stationary  frame  assunlng  U * 0,  the  following  relations  between  frames 

of  references  could  be  used 

V = U - U$  and  v0  = ' Us 

(19) 

Substituting  above  In  (16),  (17)  and  (18)  the  corresponding  stationary 
conservation  equations  reduce  to 

»o"s  ■ • U) 

(20) 

P - p0  ■ o0llsl) 

(21) 

E • (”  ^ • * (P  * ~ V) 

(22) 

Assuming  PQ  = EQ  ■ 0,  and  substituting  E from  above  Into  (*),  we  obtain: 

P = (Aj.  + Bu2  + Cy3)/(1  - ur/2)  (23) 

Similarly,  substituting  E from  (22)  Into  (15): 

t « ppr(r  + 1)(V0  - V)/2  + A + A'u  + B'u2  + C'u3  (24) 
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Finally,  substituting  P from  (23)  into  (24)  the  final  form  is  obtained  as 
Kf  = L(r  * -Bu— ^ + A + A'vj  + B'p2  + C'n3  * (25) 

where 

Ke  = the  loading  bulk  modulus 

r ■ the  Brunei  sen  coefficient 
A’  - A(r  +1)  + 2B 
B'  * B(r  + 2)  + 3C 
C'  » C(r  + 3) 

v = VO  - rv) 

cv  = ^VQ  - v)/Vq,  volume  strain  taken  positive  in  compression 

VQ  - l/pQ  = specific  volume  at  normal  conditions. 

For  the  particular  case  when  * 0 which  Implies  u * 0 and  V * VQ  it  can 

be  easily  seen  from  Equation  (25)  that  *o  = A.  Also  in  the  Gruneisen 

equation  of  state,  at  r.y  = 0,  we  take  both  the  pressure  and  internal  energy 

to  be  zero.  The  values  for  the  material  constants  of  the  explosives  used 
are  shown  in  Table  1 . 


TABLE  1 . MATERIAL  CONSTANTS  FOR  EXPLOSIVES  AND  SOIL 


Type 

pn  r A 

(9/cm3)  <Gpa) 

B 

(Gpa) 

C v 

(Gpa) 

TNT 

1.614  .737  10.367 

9.101 

138.33  .3 

TETRI L 

1.70  1.6  10.498 

17.8 

20.6  .3 

WET  TUFF 

2.0  1.5  21.77 

32.5 

18.33  — 

Since  data  to  relate  the  unloading  bulk  modulus  to  the  volumetric  strain 
were  unavailable,  the  bulk  moduli  for  unloading  were  assumed  to  be  identical 
to  the  moduli  for  loading  for  all  explosives.  The  loading  shear  modulus,  Gt, 

was  obtained  from  the  loading  bulk  modulus,  by  use  of  the  relationship, 

3<t(l  - 2v) 

G*  ~ "W 

(26) 

*Here  we  define  specific  volume,  V - I/p  or  pV 

= 1. 
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Figures  5,  6 show  the  dependence  of  bulk  and  shear  moduli  of  the  two  explosives 
represented  by  Equations  (2 5)  and  (26)  on  volumetric  strain.  Table  2 gives  the 
values  of  the  two  moduli  as  they  were  used  in  the  ADINA.  program.  ADINA  uses 
linear  interpolation  between  discrete  points. 

The  tensile  volumetric  strain  at  failure  of  -0.1  per  cent  as  given  in 
Reference  7 was  used  in  calculations  for  all  explosives  in  this  investigation. 

The  method  of  application  of  this  failure  criterion  in  the  ADINA  code  is 
throuqh  the  technique  of  superimposing  on  the  load-induced  strains,  a localized 
gravitational  pressure  sufficient  to  cause  a hydrostatic  compression  equal  in 
magnitude  to  the  tensile  failure.  When  the  total  strain  becomes  tensile  or 
negative  (as  per  convention  used  in  the  material  model  in  the  code),  a tension 
cut-off  plane  is  assumed  to  form  normal  to  the  principal  strain  direction.  The 
normal  and  shear  stiffnesses  across  this  plane  are  reduced  by  a factor  determined 
by  an  input  value.  One  or  two  additional  planes  orthogonal  to  existing  cut-off 
plane  (s)  are  allowed  to  form  if  the  strain  criterion  is  met.  The  planes  be- 
come inactive  if  compression  again  develops  in  the  normal  direction  to  the 
planes . 

The  pseudo-hydrostatic  prestrain  is  applied  by  positioning  the  vertical 
Z-coordinate  at  an  appropriate  negative  value.  The  hydrostatic  pressure 
applied  at  on  element  integration  point  is  given  for  an  element,  j,  by 


N 


V- 


(27) 


i-  i 


where 


g 

is 

the 

pe 

is 

the 

h.  . 

’J 

is 

the 

ZU 

is 

the 

N 

is 

the 

is  the  shape  function  for  node  i of  element  j. 


is  the  total  number  of  nodes  in  the  element. 


The  depth  of  the  overburden  in  terms  of  the  system  vertical  coordinate 
position  can  be  obtained  from  the  equation. 


IC  1 

0 v 


ave  gP 


(28) 


where 
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Figure  6.  Bulk  and  Shear  Moduli  vs  Volume  Strain  for 
Tetryl  Fuze  Well  Charge  of  the  TM-46  Mine 
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TABLE  2.  ADINA  INPUT  VALUES  FOR  BULK  AND 
SHEAR  MODULI  FOR  FILLER  MATERIALS 


TNT  EXPLOSIVE 


3 I 


Point  No. 

£V 

KU 

6* 

1 1 
i 

(*> 

(GPa) 

(GPa) 

(GPa) 

| 

• 

1 

0 

21 .72 

21.72 

10.62 

| 

B 

2 

1.0 

23.03 

23.03 

11.24 

! 

q 

3 

3.0 

25.65 

25.65 

12.55 
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is  the  initial  bulk  loading  modulus, 

r.y  is  the  volumetric  failure  strain,  negative  in  tension,  and 

Z is  the  negative  of  the  distance  from  the  ground  surface  to  the 

mid-plane  of  the  mine. 

C.  Soil  Simulation.  For  the  structural  response  calculations  of 
the  shallow  buried  mine,  onTy  the  top  of  the  mine  was  exposed  to  blast  pressure 
while  the  remainder  was  assumed  to  be  embedded  in  soil.  In  the  M- 1 5 mine  in 
Reference  1 an  implicit  modeling  technique  was  employed  whereby  nodal  tie 
elements  were  used  to  model  the  base  support. 

However  for  the  TM-46  mine  it  is  proposed  to  use  an  explicit  modeling 
technique  whereby  initially  two  compressible  layers  of  soil  surrounding  the 
mine  body  could  be  included  through  finite  element  discretization.  Although 
separation  of  the  mine  from  the  soil  medium  subsequent  to  the  initial  response 
and  sliding  phenomenon  could  not  be  accounted  for,  the  technique  would  be  a 
considerable  improvement  over  previous  methods  due  to  realistic  simulation 
of  the  response  which  Includes  blast  attenuation  effects.  No  simulation  of 
the  soil  was  necessary  for  the  rigid  support  calculations. 

Due  to  the  large  variety  of  soils  in  which  mines  would  be  emplaced,  it 
is  possible  only  to  select  a soil  simulation  model  which  would  be  represen- 
tative of  some  subclass  of  soils.  Thus,  a typical  shock  Hugoniot  curve  for 
wet  tuff  from  Reference  8,  9 was  selected  to  define  the  soil  element  properties. 
The  data  were  reduced  to  a high  pressure  equation  of  state  and  subsequently  to 
a dependent  formulation  of  the  bulk  moduli  on  volumetric  strain  in  order  to  be 
compatible  with  the  requirements  for  the  tension  cut-off  curve  description 
model  in  the  ADINA  code.  For  a fully  buried  mine  the  soil  modeling  could  be 
extended  in  the  region  above  the  top  of  the  mine  casing  and  only  the  top  layer 
of  the  soil  medium  could  be  pressurized  by  the  blast  load.  The  detailed  ADINA 
input  values  for  the  soil  are  shown  in  Table  2. 

P.  Simulation  of  Void  in  TM-46  Mine.  The  TM-46  mine  has  a cavity 
between  the  top  pressure  plate  and  the'  middle'  plate  covering  the  primary 
charge.  This  cavity  is  filled  with  air  which  transfers  some  load  to  the 
middle  plate  as  the  volume  of  the  cavity  is  decreased  sufficiently.  Further 
uncrimping  of  the  crimped  joint  connecting  the  side,  top  and  intermediate 
walls,  due  to  air  compression  and  angular  deflection  of  the  casing  near  the 
joint  could  conceivably  occur  resulting  in  loss  of  some  primary  explosive 
charge  and  consequent  deactivation  of  the  mine.  Although  it  is  difficult  to 
predict  uncrimping  using  the  finite  element  method,  the  code  could  be  used 
profitably  to  yield  moments  generated  and  angular  deflection  of  the  casing 
at  critical  sections.  A separate  analysis  based  on  classical  theory  could 
then  be  applied  to  predict  occurrence  of  uncrimping. 

The  airgap  Inside  the  mine  cavity  was  represented  in  the  finite  element 
model  as  a set  of  2-D  axisymmetric  fluid  elements  composed  of  an  Inviscid 
linear  compressible  material.  A constant  bulk  modulus  was  used  in  lieu  of  a 
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pressure  dependent  bulk  modulus  due  to  a lack  of  available  data  for  air. 

However,  the  primary  difficulty  with  this  model  was  that  there  was  nothing 
in  the  model  to  prevent  the  upper  plate  from  penetrating  the  middle  plate 
as  the  deformation  progressed. 

Since  the  air  was  judged  to  apply  only  a minimal  restraint  on  the  motion 
of  the  upper  plate  and  due  to  the  need  to  prevent  the  two  plates  passing 
through  one  another,  a different  model  has  been  adopted  to  simulate  contact 
and  avoid  interpenetration. 

E.  Simulation  of  Contact  in  TM-46  Mine.  Due  to  lack  of  contact 
elements  Tn  the  code  along  the  interface  between  the  top  and  intermediate 
plates  substantial  interpenetration  occurs  without  any  transfer  of  loading. 

As  a result  a major  part  of  deformation  is  confined  to  the  upper  plate, 
particularly  in  the  stair-cased  region,  which  is  clearly  unrealistic. 

Pending  modification  of  the  code  to  Include  contact  capability,  other 
alternatives  were  considered  to  overcome  the  problem  without  significantly 
altering  the  response  behavior  in  an  unrealistic  manner.  A few  alternatives 
were  eliminated  due,  either  to  a lack  of  initial  stress  Input  capability,  or 
Inability  to  vary  the  contact  element  stiffness  as  a step  function  with  axial 
compression  at  or  near  the  time  of  contact  to  prevent  overflow  at  the  interface. 
The  model  finally  selected  consists  of  fictitious  axial  truss  elements  connecting 
the  two  circular  plates.  The  material  model  for  the  trusses  is  nonlinear  and 
develops  only  a small  force  up  until  the  axial  strain  In  trusses  approaches  -1. 

At  this  strain,  a large  stiffness  is  specified  to  simulate  contact  between  the 
two  plates.  Constraints  are  applied  to  the  upper  end  of  the  trusses  to  insure 
that  its  radial  coordinate  Is  the  same  as  the  radial  coordinate  at  the  lower 
end.  Also,  the  axial  coordinate  of  upper  end  is  constrained  to  translate  with 
the  upper  plate.  These  constraints  are  necessary  to  prevent  element  rotation. 

The  model  allows  movement  only  in  the  axial  direction  and  relaxation  upon 
Initial  contact  but  it  excludes  sliding  surface  capability  which  could  conceiv- 
ably be  introduced  through  transverse  or  radial  trusses  but  would  make  the 
model  unnecessarily  complicated.  However,  for  Initial  simulation  only  areas 
with  high  probability  of  contact  have  been  considered  for  contact  simulation. 

4.  FINITE  ELEMENT  MODEL  DESCRIPTION  AND  CALCULATIONS. 

A.  Mesh  Generation.  The  finite  element  mesh  for  the  mine  was 
generated  with  the  aid  of  the  GEN3D  mesh  generator  program.  The  mine 
was  modeled  as  an  assembly  of  axlsyinnetrlc  2-D  structures  using  the  ADINA 
finite  element  code.  A six  node  QUAD  element  with  quadratic  displacement 
Interpolation  functions  In  the  direction  parallel  to  the  surface  was  used 
for  the  steel  casing.  This  type  of  element  models  the  bending  of  the  thin 
metal  casing  better  than  a four  node  QUAD.  The  explosive  components  were 
modeled  with  four-node  QUAD  elements  except  where  they  Interfaced  the  steel 
jacket.  In  which  case  a mid-side  node  was  Included  on  the  Interface  edge. 

2><2  Gaussian  Integration  points  were  used  at  each  element  for  computational 
purpose.  A total  of  304  nodes  and  157  elements  were  used  to  represent  the 
axi symmetric  model. 
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In  ADINA,  each  material  having  a distinct  material  formulation  must 
be  modeled  as  a separate  element  group.  In  general  four  major  groups  of 
elements  were  represented:  (1)  nonlinear  2-D  elements  for  the  steel  case, 

(21  nonlinear  curve  description  2-D  elements  for  the  primary  charge, 

(3)  nonlinear  2-D  elements  for  the  fuzewell  booster  charge,  and  (4)  non- 
linear truss  elements.  For  the  steel  case,  four  material  subtypes  were  used 
to  model  the  steel  properties  In  different  regions  of  the  Inhomogeneous  mine 
body.  The  material  model  for  the  casing  was  a bilinear,  von  Mises  yield 
condition,  kinematic  hardening,  2-D  ax 1 symmetric  element  model. 

B.  Time  Step  Solution.  In  ADINA,  one  has  the  choice  of  marching 
the  dynamic  solution  forward  in  time  via  explicit  or  Implicit  finite-difference 
techniques.  The  Implicit  schemes  are  unconditionally  stable  and  can  tolerate 
larger  time  step  size  resulting  In  reduction  of  computational  times.  However, 
equilibrium  iteration  and  stiffness  reformation  at  regular  intervals  are 
necessary  to  obtain  meaningful  results.  In  general,  it  is  difficult  to  make 
absolute  statements  as  to  which  Is  best  for  a given  application.  For  highly 
transient  shock  loads  such  as  shown  In  Figure  2,  it  has  been  our  experience 
that  the  explicit  method  gives  the  higher  quality  solution  for  a given 
amount  of  central  processor  computer  time.  The  subject  problem  was  run  for  a 
relatively  large  nunber  of  cycles  using  both  explicit  and  implicit  (with 
equilibrium  iterations  Included)  time  Integration  solutions.  After  a selected 
amount  of  problem  solution  time,  the  results  were  compared  for  solution 
quality.  The  explicit  solution  appeared  to  have  a smoother  variation  In  both 
displacements  and  stresses.  For  this  reason,  we  selected  the  explicit 
solution  method. 

The  time  step  used  for  the  calculations  was  determined  from  the  Courant 
stability  condition 


. atcrlt  , aLnin 
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A^m1n 
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Eymax 
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p 
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Is 

the  minimum  Courant  stability  step  size. 


the  distance  between  the  two  closest  nodes  In  the  system, 
the  Young's  modulus  for  the  stiffest  material, 
the  density  of  the  material, 

the  nunfcer  of  time  steps  which  we  wish  to  represent  the 
shock  wave  In  passing  through  the  distance  AL. 


The  value  of  Atcr^t  was  approximately  2C0  nanoseconds  and  a value  of  n of 

four  was  used,  so  that  the  time  step  for  the  central  difference  explicit  time 
Integration  method  was  50  nanoseconds. 
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C.  TM-46  Mine  Calculations . The  ADIHA  calculations  for  the  TM-46 
mine  are  in  progress.  However,  some  of  the  salient  features  of  the  model 
have  been  developed  from  progress  made  in  studies  of  the  mine  thus  far.  A 
drawing  of  the  current  mesh  configuration  is  shown  in  Figure  7. 


Figure  7.  Finite  Element  Mesh  for  the  TM-46  Mine 

Since  experimental  data  indicated  significantly  varying  material 
properties  In  different  regions  of  the  outer  steel  jacket  of  the  TM-46  mine 
several  different  sets  of  materials  properties  were  used  to  model  the  various 
steel  components  of  the  mine. 

At  the  outset  two  particular  difficulties  were  expected  to  be  encountered 
In  modeling  the  TM-46  mine.  First,  the  difference  in  stiffness  between  the 
steel  plates  and  the  air  filled  region  leads  to  numerical  problems.  The 
collapse  of  the  air  filled  region  leads  to  vhe  impact  of  the  upper  plate  on 
the  middle  plate.  This  phenomenon  needs  to  be  modeled  rather  carefully. 
Secondly,  the  thin  stepped  region  of  the  top  cover  shown  In  the  upper  right 
part  of  Figure  7 leads  to  a very  Inefficient  mechanism  for  load  transfer  from 
tne  top  cover  to  the  main  mine  body.  On  the  other  hand,  failure  of  the  top 
cover  may  not  Indicate  deactivation  of  the  mine  and  a viable^ failure, mechanism 
must  Inevitably  Involve  a failure  of  the  main  mine  body. 
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Since  the  ADINA  code  does  not  currently  have  a contact  element  to  sense 
when  the  top  cover  plate  Impacts  the  intermediate  plate,  an  approach  described 
in  Section  3E  involving  nonlinear  truss  elements  has  been  used  to  approximate 
the  Interaction  of  the  two  plates. 

Eigenfrequencies  and  mode  shapes  were  also  obtained  for  the  TM-46  mine 
in  order  to  ascertain  if  the  model  has  been  formulated  correctly  and  also 
as  an  aid  in  estimating  an  appropriate  time  step  for  the  explicit  integration 
scheme.  The  eigenfrequencies  and  associated  periods  for  the  lower  modes  are 
given  In  Table  3. 

TABLE  3.  EIGENFREQUENCIES  AND  PERIODS  FOR  THE  TM-46  MINE 


Frequency 

Period 

(cps) 

(sec) 

3041 

3.288  x 10‘4 

10466 

9.555  x 10"s 

17060 

5.859  x 10'5 

31071 

3.218  x IQ"5 

All  calculations  described  herein  used  the  total  Lagrangian  formulation 
with  a lumped  mass  matrix  with  the  exception  of  the  nodal  tie  and  truss  elements. 
The  formulations  used  for  these  were  material  nonlinearity  only  and  updated 
Lagrangian  analysis  procedure,  respectively. 

5.  DYNAMIC  RESPONSE  PREDICTIONS.  Several  modifications  to  the  ADINA 
program  were  made  to  assist  us  In  Interpreting  the  response  predictions.  These 
are  described  fully  In  Reference  1.  Due  to  the  very  large  amount  of  data 
available  from  the  ADINA  results,  search  routines  were  Incorporated  to  identify 
the  extreme  (maximum/minimum)  stresses  and  strains  In  the  steel  components. 
Additionally  for  stresses  in  the  plastic  regime  a quadratic  correction  was 
applied  to  bring  the  stress  tensor  back  to  the  von  Mlses  yield  surface  and 
allow  calculations  to  proceed  without  premature  Interruptions.  The  details  of 
these  modifications  are  also  described  In  Reference  10. 

From  the  complex  structural  construction  of  the  mine,  particularly  near 
the  stair-cased  region  in  the  top  plate  which  appeared  to  be  susceptible  to 
blast  Induced  damage  and  failure,  it  was  decided  to  model  this  region  very 
accurately  using  shadowgraph  measurements.  The  elements  and  nodes  In  the  top 
plate  were  configured  to  line  up  In  the  axial  direction  with  those  In  the 
intermediate  plate  to  facilitate  inclusion  of  contact  elements  and  minimize 
element  rotation. 

The  Initial  calculations  for  the  mine  on  a rigid  base  support  were  run 
.^without  the  contact  elements  to  approximately  4000  cycles  corresponding  to  an 
elapsed  time  of  0.2  milliseconds.  The  purpose  of  the  initial  runs  were  to 
determine  the  regions  and  the  time  of  Initial  contact  between  the  two  plat’s. 

The  results  indicated  occurrence  of  Initial  contact  at  .06  milliseconds  near 
the  crimped  region  of  the  mine  body  and  the  stair  cased  region  of  the  top  cover. 
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Failures  were  predicted  in  several  elements  in  the  stepped  regions 
accompanied  by  considerable  plastic  flow.  Stresses  were  significantly  high 
particularly  near  sharp  corners  due  to  acute  stress  concentration  effects. 

A major  part  of  the  response  was  evident  in  the  top  plate  while  the  rest  of 
the  mine  body  did  not  exhibit  any  appreciable  deformation.  The  deformation 
of  the  mine  confined  mainly  in  the  area  of  the  top  cover  plate  may  not  be 
very  realistic  since  the  mine  supported  on  a rigid  roller  base  has  unrestrained 
sidewalls  where  appreciable  deformation  is  expected  but  unrealized  due  to 
inefficient  load  transfer  mechanism  between  the  top  and  intermediate  walls.  One 
of  the  chief  difficulties  encountered  has  been  in  trying  to  provide  the  appro- 
priate model  for  the  interaction  of  the  top  cover  plate  on  the  middle  plate.  We 
have  used  the  improvisation  of  nonlinear  truss  elements  (see  Figure  7)  to 
simulate  the  impact  of  these  two  components.  Since  a sudden  stiffening  of  the 
trusses  near  the  time  of  contact  generated  spurious  transients  and  was  in- 
sufficient to  arrest  penetration  of  the  intermediate  plate  due  to  inertia 
effects,  earlier  stiffening  of  the  trusses  in  a gradual  manner  initially  and 
at  an  accelerated  rate  subsequently  as  shown  in  Figure  8 was  resorted  to.  The 
trusses  developed  significantly  high  stresses  and  internal  forces  in  the  region 
where  contact  first  occurred.  Subsequently  upon  impact  the  separation  of  plates 
caused  the  stresses  in  the  affected  trusses  to  be  released  while  the  next  set 
of  trusses  in  the  neighboring  region  approached  impact  conditions  and  developed 
high  stresses  which  were  then  released  as  plates  separated  and  the  traveling 
wave  propagated  radially  inward  until  the  initial  contact  process  wns  completed. 
Although  the  response  of  the  trusses  appeared  to  be  realistic,  some  overflow 
occurred  in  the  corner  region  and  further  stiffening  was  deemed  to  be  necessary 
to  prevent,  interpenetration  and  obtain  meaningful  results. 

A typical  response  of  the  system  at  an  early  time  is  shown  in  Figure  9. 

In  this  figure  the  dotted  lines  represent  the  undeformed  or  original  configu- 
ration before  imposition  of  the  blast  load.  The  vertical  lines  between  the 
top  cover  plate  and  the  middle  plate  represent  the  nonlinear  truss  elements. 
Currently,  the  calculation  has  not  proceeded  to  the  point  where  any  failure 
of  the  main  mine  body  can  occur.  However,  failure  of  some  casing  elements  In 
the  stepped  region  has  been  observed  at  several  time  periods  beginning  at 
.02  ms  and  continuing  beyond  ,1  ms  when  formation  of  a zone  of  rupture  was 
indicated  In  the  top  cover.  However,  the  main  mine  body  did  not  show  any 
evidence  of  failure  due  to  the  inefficient  load  transfer  mechanism  between 
the  top  and  intermediate  plates.  The  model  of  this  mine  is  still  evolving 
and  Is  undergoing  modifications  to  improve  accuracy  and  reliability  of 
prediction. 

6.  CONCLUSIONS.  The  soviet  TM-46  mine  because  of  its  complex  stair- 
cased,  double-walle<f  construction  was  somewhat  difficult  to  model  accurately, 
to  analyse  and  to  deactivate.  Failure  of  a particular  region  of  the  mine 
casing  itself  might  not  be  sufficient  to  indicate  overall  failure  of  the  mine. 

A conservative  approach  was  therefore  taken  to  model  the  mine  by  eliminating 
nonessential  details  and  by  requiring  defeat  of  the  main  mine  body  through 
case  rupture. 

The  explicit  time  integration  method  appeared  to  be  advantageous  for  the 
shock  loaded  mines  due  to  smoothness  of  stresses  and  strains  as  a function  of 
time.  However,  second  order  corrections  were  necessary  to  assure  that  the  stress 
state  remained  on  the  yield  surface  during  plastic  flow. 
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Figure  8.  Input  Values  for  Stiffness  Variation  of  Trus 
Fleaents  Used  in  Hie  ADINA  Code 


The  parts  of  the  outer  steel  jacket  of  the  TM-46  mine  which  are  work 
hardened  in  the  deep  drawing  metal  forming  operation  have  significantly  vary- 
ing materials  properties.  These  variations  in  stress-strain  relations  must 
be  measured  and  modeled  carefully  since  they  directly  affect  mine  failure 
under  blast  loads. 

The  contact  problem  arising  from  interpenetration  of  plates  results  in 
inefficient  load  transfer  mechanism  between  the  cover  and  the  intermediate 
plate.  Truss  contact  elements  may  be  used  profitably  to  overcome  this  short- 
coming provided  nodal  constraints  are  used  to  avoid  truss  rotations. 

Linear  fluid  elements  can  be  used  successfully  to  represent  interaction 
of  air  in  voids  with  structural  response.  However,  in  TM-46  mine  it's  effect 
on  frequencies  and  mode  shapes  is  very  minimal. 

The  initial  deformation  of  the  TM-46  mine  was  limited  to  the  response  of 
the  top  cover  plate.  The  finite  element  modelinq  with  the  ADINA  program  has 
presented  some  difficulty  in  describing  accurately  the  impact  of  the  top  plate 
on  the  intermediate  plate  which  contains  the  primary  TNT  charqe.  From  past 
experience  main  body  failure  is  expected  to  occur  at  later  times  near  the 
corners  of  the  fuzewell  due  to  a change  of  thickness  and  a sharp  radius 
resulting  in  significantly  high  stress  concentration  and  case  rupture.  This 
study  is  still  in  progress. 
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